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XMM-Newton observations of the merging galaxy cluster 
CIZA J2242. 8+5301 
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ABSTRACT 

We studied the intracluster medium of the galaxy cluster CIZA J2242. 8+5301 using 
deep XMM-Newton observations. The cluster hosts a remarkable 2-Mpc long, ~ 50-kpc 
wide radio relic that has been nicknamed the "Sausage" . A smaller, more irregular 
counter-relic is also present, along with a faint giant radio halo. We analysed the 
distribution of the ICM physical properties, and searched for shocks by trying to 
identify density and temperature discontinuities. East of the southern relic, we find 
evidence of shock compression corresponding to a Mach number of ~ 1.3, and speculate 
that the shock extends beyond the length of the radio structure. The ICM temperature 
increases at the northern relic. More puzzling, we find a "wall" of hot gas east of the 
cluster centre. A partial elliptical ring of hot plasma appears to be present around the 
merger. While radio observations and numerical simulations predict a simple merger 
geometry, the X-ray results point towards a more complex merger scenario. 
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1 INTRODUCTION 

Galaxy clusters form hierarchically, via mergers with other 
clusters and groups of galaxies, and also by accretion of 
gas from the intergalactic medium (IGM). Mergers may 
trigger shocks and turbulence in the intracluster medium 
(ICM), and extended sources of diffuse synchrotron emission 
are typically observed in clusters that show evidence for a 
merger (e.g. iBuotd l200ll: iGovoni et alj 12004 ; IVenturi et all 
12008 ; ICassano et al.ll2oTol : Ivan Weeren et al.ll201lh . 



It has been proposed that relics trace shock waves in- 
jected into the ICM during a major merger. The commonly 
accepted theory is that in the presence of a magnetic field, 
particles are (re-)accelerated at the shock to relativistic ener- 
gies and emit synchrotron radiation. This emission is visible 
in the radio as arc-shaped relics at the cluster periphery. 

Radio observations reveal the position of the relic, while 
X-ray observations can locate the shock wave by identifying 
surface brightness discontinuities, or temperature, pressure, 
and entropy changes within the ICM. Therefore, multiwave- 
length observations of merging galaxy clusters harboring 
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relics are necessary in order to test the radio relic-shock 
wave hypothesis. Moreover, such observations provide infor- 
mation about the merger dynamics, the processes behind 
relic formation, and the intracluster magnetic fields in the 
region of the relic. At the same time, they allow us to re- 
fine models of shock acceleration in merging clusters, where 
the low Mach numbers suggest an efficiency inadequate to 
explain the existence of radi o relics by direct shock acceler- 
ation (e.g. jKang et alj|2007h . 

So far, only a handful of shock fronts have been 
found that exhibit, both, a sharp gas density edge 
and an unambiguous te mperature jump: in th e "Bul- 
let Cluster", 1E065756 jMarkevitch et~al] |2002h. Abell 
520 (|Markevitch et al.l 120051). Abell 214 6 (|Russell et ail 
l2010h. Abell 3667 dFinoguenov et al.ll2010h . and Abell 754 
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Of these, only the shocks in Abell 
3667 and Abell 754 are associated with radio relics, even 
though there are about 50 known relics. Such discoveries 
are rare because the merger has to be observed at the time 
when the shock has not yet moved to the low-brightness out- 
skirts and is propagating nearly in the plane of the sky in 
order to give a clear view of the shock discontinuity. 

Here, we present a combined X-ray and radio analy- 
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sis of the merging galaxy cluster CIZA J2242. 8+5301 (z = 
0.1921), an X-ray luminous cluster (Lx = 6.8 x 10 44 erg/s in 
the ROSAT band 0.1 - 2.4 keV) which hosts a double radio 
relic, and an extended, faint radio halo. The properties of 
the re lics have been analysed in detail by Ivan Weeren et al.l 
(|20ld ); in this introduction, we present a summary of their 
findings. The relics are oriented along the N-S direction. 
The northern relic is the most spectacular, with a length 
of 2 Mpc, and a width of only 55 kpc, narrower than any 
other known relic. At this relic, deep radio observations have 
revealed a spectral index gradient towards the centre, with 
the spectral index changing from —0.6 to —2.0 (T v cx ^~ Q , 
where T v is the radio flux at frequency v, and a is the spec- 
tral index). Furthermore, the relic is strongly polarized, at 
a level of ~ 50%, wi th magnetic field vector s aligned with 
the relic (see fig. 3 in Ivan Weeren et all (|2010t )). These radio 
properties provide clear evidence for shock acceleration and 
spe ctral ageing at an outwar d moving shock. 

Ivan Weeren et aT] f|201 lh used hydrodynamical simula- 
tions of binary cluster mergers to constrain the merger ge- 
ometry in CIZA J2242. 8+5301, based on the observed radio 
structures. They found that the cluster is most likely under- 
going a merger in the plane of the sky, between two clusters 
of nearly equal mass (mass ratio of about 2:1) that collided 
with a small impact parameter (< 400 kpc). In the following, 
we explore the spatial and spectral properties of the ICM us- 
ing a deep XMM-Newton observation of the cluster, and we 
compare the predicted merger scenario with that implied by 
the X-ray results, in order to further refine our understand- 
ing of the merger geometry that led to the formation of the 
spectacular double-relic system. 

The paper is organized as follows: Section[2]presents the 
observations and the data reduction. Section [3] details our 
analysis, while in Section U we analyse the effect of system- 
atic uncertainties on the spectral measurements. Sections [5] 
[U and [7] discuss the results. A summary of our findings is 
given in Section [5] 

We assume a flat ACDM universe with Ho = 70 
kms^Mpc -1 , O m = 0.3, and fi A = 0.7. At the redshift 
of the cluster, 1 arcmin corresponds to 192 kpc. Through- 
out the paper, quoted errors are la statistical errors, unless 
stated otherwise. Normalizations of all the spectral compo- 
nents are given in default Xspec units (i.e., cm~ J for APEC 
components, and photons keV" 1 cm -2 s _1 , measured at 1 
keV, for power-law components) per square arcmin. Tem- 
peratures are in units of keV. To enhance the clarity of the 
spectral plots, all spectra have been rebinned in Xspec to a 
minimum significance of 3cr, while setting a limit of 30 for 
the maximum number of combined adjancent bins. 
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Figure 1. Exposure-corrected and instrumental background- 
subtracted XMM-Newton EPIC surface brightness map of the 
cluster, in the energy band 0.5 — 4 keV. The image was binned 
by 2, and adaptively smoothed with a minimum of 100 counts 
per bin. Overlaid are WSRT 1.4 GHz radio contours, drawn at 
[1, 5, 10, ...] X 100 /ijy/beam. 



with 64 channels per sub-band. The data were calibrated 
using the CASAQ package. 

The reduction consisted of standard gain and bandpass 
calibration using a flux calibrator observed at the start and 
end of a run. Time ranges for antennas affected by shad- 
owing were removed and some data were flagged due to 
radio frequency interference. The fluxes for the cal i brator 
sources were set according to the Perlev fc Tavlorl (| 19991 ) 



extension to the iBaars et al. (1977) scale. Several rounds 
of phase self-calibration, followed by two rounds of ampli- 
tude and phase self-cali bration were carried out. An image 
with robust weighting |Briggsl ll995T ) of was made. The 
image was cleaned using manually placed clean boxes. The 
resolution of this image is 16.5" x 12.9" and the noise is 
19 /iJy beam -1 . 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Radio 

CIZA J2242.8+5301 was observed with the Westerbork Syn- 
thesis Radio Telescope (WSRT) in the 21 cm band. The ob- 
servations were spread out between March and November 
2009. The total time on target time was about 30 hrs. A 
bandwidth 160 MHz of was used, evenly divided over 8 sub- 
bands. All four linear polarization products were recored 



2.2 X-ray 

CIZA J2242. 8+5301 was observed on December 13-15, 2010 
for 127 ks with the three EPIC cameras of XMM-Newton. 
The observations were performed in full-frame mode, using 
the medium filter. 

Data reduction for all cameras used the Extended 
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Table 1. Residual soft proton contamination for the MOS and 
pn event files, based on the 7£sp parameter. We list the val- 
ues calculated in the hard band , i.e. 10 — 12 keV fo r MOS , and 
12 - 14 keV for pn Follo wing |Pe Luca fc Molendil j2004) and 
iLeccardi fc Molendil fcOOSft , we also calculated 7^.gp in the energy 
bands 8 — 12 and 6 — 12 keV, respectively, for the MOS event files. 



Energy band 


MOS 1 


MOS 2 


pn 


6-12 keV 


1.10 ±0.024 


1.07 ±0.022 




8-12 keV 


1.09 ±0.028 


1.10 ±0.026 




10 - 12 keV 


1.10 ±0.037 


1.08 ±0.033 




12 - 14 keV 






1.10 ±0.044 



Source Analysis Software (esas) integrated in the XMM- 
Newton Science Analysis System (sas) version 12.0.1. Raw 
event files were created from the Observation Data Files 
(ODF) using the xmm-esas routines emchain and epchain. 
As XMM-Newton orbits the Earth, protons with energies 
less than a few hundred keV can interact with the instru- 
ments, and be funneled towards the detectors. These soft 
protons (SPs) will create a time- variable (flaring) instrumen- 
tal background component inside the field of view (FOV). 
The data was filtered of SP flares using the mos-filter 
and pn-filter routines, which apply a 1.5a cutoff around 
the center of a Gaussian fitted to the EPIC lightcurve his- 
tograms. We then checked each filtered event file for residual 
soft proton contamination. This was done by comparing the 
EPIC count rates inside the field of view (thus, affected by 
possible residual soft protons) to the corresponding count 
rates in the unexposed corners of the detectors. The count 
rates were calculated in the energy bands 10—12 keV (MOS) 
and 12 — 14 keV (pn), in which the detectors are essen- 
tially insensitive to X-ray photons. The higher the count 
rate ratio inside-to-outside the FOV (IZsp) in these energy 
bands, the higher the lev el of residual SP contamination. 
|Pe Luca fc Molendil (|2004l ) found that count rate ratios be- 
low 1.15 indicate event files that are hardly contaminated by 
residual soft protons. According to the same authors, count 
rate ratios between 1.15 and 1.3 correspond to event files 
slightly contaminated by soft protons, while ratios above 1.3 
signify that the event files are significantly contaminated by 
flares. The values of T^sp following the mos-filter run are 
shown in Table[T] For the pn data, the flare-filtered event file 
created by the pn-filter routine had "Rsp = 1.27 ± 0.049. 
Therefore, we created pn light curves and count rate his- 
tograms in the energy bands 2 — 5 and 12 — 14 keV (most 
sensitive to flares), and gradually lowered the count rate 
threshold of good-time intervals down to the maximum level 
at which 7^sp was below 1.15. After filtering, the total good- 
time intervals were 60 ks for MOS1, 62 ks for MOS2, and 25 
ks for pn. Table [T] lists the 7?.sp values for the three filtered 
EPIC event files. 

Contaminating point-like sources within the FOV were 
identified with the xmm-esas routine cheese, using a flux 
threshold of 3 x 10~ 15 erg cm -2 s _1 . The output files were 
then checked for spurious detections and undetected sources. 
All point-like sources were excluded from the analysis. 

Spectra and images were extracted using the rou- 
tines mos-spectra and pn-spectra. The resulting data, 



along with the xmm-esas CalDB filefl describing the quies- 
cent particle background, were then used by mos_back and 
pn_back to create images and spectra of the quiescent parti- 
cle background (QPB). We used the XMM-Newton CalDB 
files released in June 2012. 

The pn data was corrected for out-of-time (OoT) events. 
OoT events are recorded during the readout time of the 
CCDs, and are therefore assigned an incorrect RAWY posi- 
tion. For the full-frame mode, the ratio of the readout time 
to the total integration time is 6.3% for the pn detector 
(only 0.35% for MOS, for which no correction is necessary). 
Therefore, the xmm-esas routines create spectra and images 
of the OoT events, scale them by 6.3%, and subtract them 
from the pn images and spectra. 

CCD #6 of MOS1 became unoperational in 2005, and 
is automatically excluded by xmm-esas. 



3 DATA ANALYSIS 

3.1 Imaging analysis 

To search for surface brightness edges, we created an image 
in the energy band 0.5 — 4 keV, excluding the energy range 

1.2 — 1.9 keV, which contains emission from the instrumental 
Al Ka and Si Ka lines (see Section l3~^l for details). At higher 
energies, the emission is dominated by background, hence 
our upper-energy cutoff. 

Individual detector images were created using the 
tasks mos-spectra and pn-spectra, and combined with the 
comb routine. We used the xmm-esas task adapt_900 to 
create an instrumental-background-subtracted, vignetting- 
corrected EPIC image, binned by 2, and adaptively 
smoothed with a minimum of 100 counts per bin. Due to the 
size of the XMM-Newton FOV, we also limited the image to 
a circle of radius 12 arcmin around each detector's centre. 
The resulting image is shown in Figures [T] and [2] The X- 
ray morphology is strongly elongated along the north-south 
direction, evidence that the cluster is undergoing a merger 
along this axis. At ~ 170 kpc from the centre there is an 
arc of X-ray emission. Furthermore, a "smudge" extending 
on a length of about 500 kpc is aligned with the relic on its 
easternmost side. To the east, the tip of the smudge coin- 
cides with the position of a n AGN, seen as a brigh t radio 
point-like source in Figure [T1 (|van Weeren et "al]|2010l ). Inter- 
estingly however, the AGN jets are oriented perpendicular 
to the relic, i.e. also perpendicular to the smudge, so it is not 
clear whether the X-ray emission originates from the AGN, 
or is instead associated with a density and temperature en- 
hancement at a shock tracing the radio relic. The arc and 
the "smudge" feature are shown in Figure [2] 

We examined the surface brightness of the cluster in 
six directions from the centre, as shown in Figure [3] where 
the centre was defined as the geometrical centre of the 
double-relic system. For this, we created an instrumental- 
background-subtracted, vignetting-corrected EPIC surface 
brightness map in the energy band 0.5 — 4 keV, binned by 2. 
To avoid large vignetting-correction uncertainties near the 
edges of the detectors, we limited the profiles to the inner 24 
arcmin (a circle centred on each detector centre, and having 

2 ftp: / / xmm.esac.esa.int / pub / ccf /constituents/extras / esas_caldb / 
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Figure 2. Top: Exposure-corrected and instrumental background-subtracted XMM-Newton EPIC surface brightness map of the cluster, 
in the energy band 0.5 — 4 keV. The image was binned by 2, and adaptively smoothed with a minimum of 100 counts per bin. Bottom: 
(A) Zoom-in on the "smudge" feature. (B) Zoom-in on the arc. The bottom two images have a side length of 0.8 Mpc. 
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a radius of 12 arcmin). We chose the northern pie sector 
in such a manner as to cover the northern relic, but avoid 
the northeastern AGN and the observed X-ray "smudge" 
(see Figure [2] and also Section [6]). East of the northern sec- 
tor, we selected a narrower one which covers the "smudge" 
region. The other four sectors' angular widths were chosen 
such that surface brightness contours within each of them 
are roughly circular far from the cluster centre. We assumed 
a constant sky background surface brightness, and calcu- 
lated its contribution from the crescent-moon shaped back- 
ground region shown in Figure [5] The background surface 
brightness was then subtracted from the map. In each sec- 
tor, the average surface brightness was calculated in partial 
annuli of constant signal-to- noise ratio (SNR), SNR = 5. 
All errors were propagated in quadrature. The background- 
subtracted surface brightness profiles are shown in Figure [4] 
We note that close to the centre, the complex substructure 
is not well-described by means of surface brightness profiles. 
However, our aim is to search for signatures of shock com- 
pression in the outskirts, where there is little substructure, 
and the merger's ellipticity can be safely ignored inside the 
selected sectors. 

Beyond the northern relic, there is a single bin of SNR ~ 
1.8 (not shown in Figure [4]) and surface brightness 0.08 ± 
0.05, so a factor of ~ 7 lower than the surface brightness 
in the (putative) post-shock region. Unfortunately, the ICM 
signal beyond the relic is too low to allow us to model a 
putative shock with a dou ble power-law densit y profile, as 
it is customarily done (e.g. iMacario" et al.ll201lf ). 

Near the southern relic, we observe a surface brightness 
jump by a factor of ~ 2 — 3. Here, the small number of 
pre-shock bins with good SNR make modelling the density 
beyond the relic difficult. Instead, we model the pre-shock 
density assuming a bet a-model descr i bed b y the average /3 
power index derived bv lEckert et~aH (|2012h for a sample of 
31 galaxy cluster observed with ROSAT. The method and 
results are presented in detail in Section [5] 

In the narrow bin emcompassing the "smudge", there 
is a clear peak in the surface brightness at approximately 
1.3 Mpc from the centre, coinciding with the position of the 
observed "smudge" . The origin of this emission is discussed 
in Section [6] 
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3.2 Background and foreground modelling 

The low surface brightness in cluster outskirts makes it es- 
sential to model accurately the different background and 
foreground (BF) components. Typically, there are four BF 
components that need to be considered: a quiescent par- 
ticle background (QPB), produced by the interaction of 
high-energy particles with the detectors and other satel- 
lite components; unresolved background sources; the Galac- 
tic Halo (GH); and the Local Hot Bubble (LHB). How- 
ever, CIZA J2242. 8+5301 is located at a low Galactic lati- 
tude, |6| « 5°, where add itional foreground c omponents can 
be present. For example, iMasui et al.1 (|2009D detected ther- 
mal foreground emission with T ~ 0.8 keV in the Galac- 
tic direction (235°, 0°), which they attribute most likely to 
a population of unresolved dM stars, w hose number falls 
quickl y with increasing Galactic latitude. ISimionescu et all 
(|201lf ) found a foreground component with T w 0.6 keV to- 



Figure 4. Background-subtracted surface brightness profiles in 
the sectors shown in Figure [3] Only bins with a SNR > 3 are 
shown. Error bars are Poisson errors propagated in quadrature. 
Dotted lines show the positions of the relics based on the WSRT 
1.4 GHz observations. Profiles extend up to 1.5 Mpc. 

wards the P e rseus cluster, at Galactic latitude |b| ~ 13°n 
lEckert et all <|201lf ) compared Suzaku (jGeorge et al.ll2009T l 
and ROSAT observations of the galaxy cluster PKS 0745- 
191 (|6| ~ 3°), and found discrepancies between the surface 
brightness profiles measured with the two satellites; they at- 
tribute these discrepancies to additional backg round compo- 
nents that were not included in the analysis of iGeorge et alj 

3 ISimionescu et all l|201ll ) calls this component the "hot fore- 
ground" (HF), which is also how we will refer to it in this work. 
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Figure 3. Sectors used for creating surface brightness profiles. Dashed circles mark distances of 0.5, 1.0, and 1.5 Mpc from the centre, 
and are shown for easier comparison with the profiles in Figure [4] 



(|2009h . We note that the P KS 0745-191 Suzak u data has 
been re-analysed recently bv lWalker et~aI1 (|2012h using new 
background observations near the cluster, and it was found 
that the inclusion of an additional background component 
with T ps 0.5 keV provides the best fit to the background and 
eliminates the discrepancy between different instruments. 

To model the background near CIZA J2242. 8+5301, as 
well as all the other spectra presented in this paper, we used 
version 12.7.1 of Xspec. EPIC spectra were extracted from 
the background region shown in Figure [5] and grouped to 
a minimum of 30 counts/bin. To better constrain the soft 



components of the BF model, we also used four ROSAT 
spectra extracted from an annulus around the cluster centre, 
and from three circles NW and SE of the cluster centre. 
We avoided the NE direction, where another galaxy cluster 
is visible in the ROSAT image, and the S-SW directions, 
where the X-ray count rate is slightly higher. The positions 
and extent of the selected regions are summarized in Table 
H while Figure [5] shows their location on the ROSAT PSPC 
contour map. 

xmm-esas uses the CalDB files to model the QPB. How- 
ever, it excises from the modelled QPB spectra the strong 
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Figure 5. ROSAT PSPC contour map, showing equidistant con- 
tours. Overlaid in solid lines are the regions used for extract- 
ing ROSAT All-Sky Survey (RASS) background spectra. The 
crescent-moon shaped region in dashed lines shows the region 
used for extracting the XMM-Newton background spectra. 



Table 2. Regions used for extracting ROSAT background spec- 
tra. Region names correspond to the labels in Figure [5] 



Region 




RA 




DEC 


Radius (deg) 


Rl 


22 h 


39 m 


23 s 


+53° 19' 27" 


0.25 


R2 


22 h 


46 m 




+52° 51' 40" 


0.25 


R3 


22 h 


42 m 


01 s 


+53° 34' 53" 


0.22 


R4 


22 h 


42 m 


45 s 


+53° 02' 40" 


0.26, 0.33 



Table 3. Best-fit background and foreground parameters. 



Component 


r 


T x 


JVx 


LHB 




0.08t 


(3.4 ±0.29) x 10~ 7 


GH 




Q -, 4 +0.0027 
U,J - -0. 0(124 


1 fi+ 095 v TO- 5 
i -°-0.094 X iU 


HF 




n fi9 + ().02:i 
u -° -0.022 


(8.7 ±0.59) x TO" 7 


CXBtt 


1.41+ 




1 n+ 046 x W~ 6 

i ' U -0.045 X iU 



f frozen 

ff for comparison, iDe Luca fc Molendil ||2004TI obtain a total 
CXB normalization (i.e., including both resolved and unresolved 
sources) of (9.8±0.35) X 10" 7 photons keV -1 cm 2 s 1 arcmin 2 
at 1 kcV 



3.3 Temperature and pressure distribution 



AI Ka (~ 1.5 keV; present in the MOS and pn spectra), 
Si Ka (« 1.75 keV; present in the MOS spectra) and Cu 
(~ 8 keV; present in the pn spectra) fluorescent instrumen- 
tal lines. Therefore, these lines must either be added later to 
the spectral model describing cluster emission, or the energy 
sub-bands containing these lines need to be ignored in the 
fitting process. Fitting the Cu lines is difficult. Moreover, 
the pn data also has a strong low-energy tail at energies be- 
low 0.3 — 0.5 keV. In consequence, we fitted the spectra in 
the energy band 0.5 — 7.0 keV, excluding the 1.2 — 1.9 keV 
MOS data points, and the 1.2 — 1.7 keV pn data points. 

Foregroun d emission was m odelled with the sum of 
three APEC (|Foster et all l2012t ) components, describing 
emission from collisionally-ionized diffuse gas: an unab- 
sorbed component, corresponding to LHB emission, and 
two absorbed components, corresponding to GH and HF 
emission, respectively. The LHB temperature was fixed to 
0.0 8 keV, which is the aver age of the tempera t ures found 
by ISidher et ail |l99fj ) and iKuntz fc Snowdenl |2000[) . We 
used the abundance table of lAnders fc Greve ssc ( 198S ), and 
fixed the abundances of the foreground components to so- 
lar value. The redshifts were fixed to 0. The other APEC 
spectral parameters were free, but coupled between spec- 
tra. The cosmic X-ray background (CXB), consisting of 
emission from unresolved cosmological sources, was mod- 
elled with an a bsorbed power-law with a fixed photon in- 
dex r = 1.41 (|De Luca fe Molendil l20ui ) . The CXB nor- 
malization was free, but coupled between the XMM-Newton 
spectra. For the ROSAT spectra, the CXB normalization 
was fixed to 8.85 x 10~ 7 photons keV" 1 cm~ 2 s" 1 arcmin -2 
jMoretti et al.ll2003l ), as point sources were not excluded 
from these spectra. Photoelectric absorption was modelled 
with a PHABS component, and we used the p hotoelectric 
absorption cross-sections of IVerner et all (|l99fi ). 



We use d the contour binning algorithm developed by 
ISandersI |2006l r1 to extract spectra with a total of ~ 3600 
ICM plus BF counts. The code bins the data following the 
surface brightness contours of the map shown in Figure [2] 
This approach is better suited than, for example, extract- 
ing spectra from annuli or partial annuli around the centre, 
given that the cluster has an extremely disturbed morphol- 
ogy, and that we are searching for density/temperature dis- 
continuities in the ICM. We used a smoothing signal -to-noise 
of 30, and a geometric constraint value of 1.25 (see ISandersI 
l|2006h for details). The small constraint value ensures that 
the bins do not become narrower than the XMM-Newton 
PSF. 

All spectra, including those of the background, were fit- 
ted simultaneously, in order to maximize the number of con- 
straints and correctly propagate the statistical errors. The 
XMM-Newton spectra were normalized by their respective 
geometrical areas (in units of sq arcmin) when defining the 
spectral model. The ROSAT spectra are already in units of 
per sq arcmin. It was assumed that the ICM temperature in 
each bin is well-approximated by a single-temperature spec- 
trum, which we modelled using an absorbed APEC compo- 
nent added to the background model described in the pre- 
vious section. The ICM normalization was fixed to zero for 
the background spectra. 

In the direction of the cluster (RA = 22 h 42 m 53 s , 
DEC = +53° 01' 05"), in a circle of radius 1 degree, the 
weighted average atomic hydrogen column density (Ah) 
listed in t he Leiden/Argentine /Bonn (LAB) Survey of 
Galactic Hi (|Kalberla et al.ll2005l ) is 3.22 x 10 21 cm -2 . How- 
ever, there is a large scatter in the LAB A^h values, with 
Galactic hydrogen column densities taking values in the 
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Figure 6. Temperature and pseudo-pressure maps of the cluster. Spectra were extracted from bins that follow the surface brightness 
contours, and have a total of ~ 3600 counts after the subtraction of the instrumental background. The pressure was calculated as the 
product between the temperature and the square root of the normalization of the spectral component describing ICM emission. Overlaid 
are the radio contours also shown in Figure [T] 



range (2.53 - 3.65) x 10 21 cm" 2 . Therefore, the X-ray col- 
umn density for each ICM spectrum was a free parameter 
in the fit. For the background spectra, there is some de- 
generacy between the column density and the foreground 
components, so the X-ray column density was fixed to the 
weighted average LAB value. We discuss in Section [4] how 
this affects our best-fit ICM parameters. The fit to the spec- 
tra had a x 2 /d.o.f. of 4909.13/4358. A spectral model with- 
out HF emission provided a poorer fit to the data, with 
X 2 /d.o.f. = 5104.75/4360. We also ran an F-test to check 
that the HF component is a reasonable addition to the 
model, and obtained a F-test probability ~ 10~ 37 . There- 
fore, a spectral model that includes HF emission provides 
the best fit to the data, as it was already discussed in more 
detail in Ogrean et al. (2012). 

Figure [TJ] shows the resulting temperature and pseudo- 
pressure maps. The pseudo-pressure was calculated as the 
product between the temperature and the square root of the 
normalization. The fit to the background spectrum is shown 
in Figure and in Table [3] we list the best-fit background 
parameters. Table I All of the Appendix gives the tempera- 
ture and normalization in each bin, together with the cor- 
responding lcr errors. Figure [ATI in the Appendix shows the 
bin map, while Figure lA"2l presents the fitted spectra of all 
the bins. For comparison, we also fitted the spectra using 
an X-ray column density fixed to the average LAB value. 
The results of this fit are listed in Tables IB 1 1 and IAU the 



SKY BACKGROUND 



0.01 




1 2 

Energy (keV) 

Figure 7. Fit to the EPIC background spectrum extracted from 
the region marked in Figure[5] MOS1 and MOS2 data are plotted 
in black and red, respectively, while pn is plotted in green. The 
bottom panncl shows the ratio of the data to the model. 



corresponding spectra are shown in Figure [B2l while Figure 
IB1I shows the temperature and pseudo- pressure maps. 
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Figure 8. The effect of sky background systematic errors and TVh uncertainties on the best-fit temperatures and normalizations. The 
minimum and maximum parameters for each bin were found by varying the sky background spectrum parameters within their MIN— 1<7 
and MAX+lc ranges, where the MIN— la and MAX+lcr values were the best-fit parameters and corresponding ±lcr errors obtained, 
respectively, for the minimum and maximum Nh listed in the LAB Survey. See the text for details. 



4 SYSTEMATIC UNCERTAINTIES 

4.1 TVh uncertainties and sky background 
systematic errors 

Throughout the paper, the errors quoted are la statistical 
errors. Therefore, the effects of systematic uncertainties on 
the fitted ICM spectra have so far been ignored. Further- 
more, in fitting the spectra, we assumed the X-ray column 
density at the position of all five background regions was 
equal to the weighted average Hi column density listed in 
the LAB Survey, Nh = 3.22 x fO 21 cm -2 , calculated in 
a circle of radius 1 degree centred on the cluster coordi- 
nates in the NASA/IPAC Extragalactic Database (NED), 
RA = 22 h 42 m 53 s , DEC = +53° 01' 05". Because of the 
cluster's low Galactic latitude, there is a scatter of 20 — 25% 
in the Galactic Hi column density within this l°-radius cone. 
Local variations in Nh could affect our spectral measure- 
ments. Here, we examine the effect of systematic uncertain- 
ties and A'h fluctuations on the temperature and pressure 
distribution. To do this, we refitted the background spectra 
alone, with the X-ray column density fixed to the minimum 
and maximum Galactic hydrogen column densities listed in 
the LAB Survey within 1 degree of the cluster's centre, i.e. 
Nh, mm = 2.53 x 10 21 cm" 2 and JV H , max = 3.65 x 10 21 cm" 2 . 



Using the MIN — la and MAX+lcr value of each 6 BF param- 
eters, we then constructed 12 background models, each hav- 
ing one BF parameter fixed to its MIN — la or MAX+lcr value, 
while the other parameters were frozen to the best-fit val- 
ues in Table [3] Finally, the 16 ICM spectra were refitted 
using each of the 12 background models and a free X-ray 
column density. Figure[S]presents a comparison between the 
maximum and minimum best-fit temperatures and normal- 
izations, and the best-fit values listed in Table IA1I Within 
the la statistical errors of each measurement, the results are 
consistent for all bins, with the exception of bins 13 and 14. 
However, the large variations introduced by A'h uncertain- 
ties on the best-fit parameters of bins 13 and 14 does not 
change any of the main results presented in this paper. 



In conclusion, our main spectral results are robust, in 
spite of uncertainties on the X-ray column density in the 
direction of the background regions. Even by varying the Nh 
corresponding to the background regions by more than 1 x 
10 21 cm" 2 , the temperature and normalization distributions 
are qualitatively similar to those shown in Figure [6] 
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Figure 9. The effect of a 5% instrumental background uncertainty on the best-fit temperatures and normalizations. 



4.2 Uncertainties on the QPB 

We also examined the spectral effects of a 5% uncertainty on 
the normalization of the quiescent particle background. We 
first varied the BACKSCAL keyword of the instrumental 
background spectra associated with the sky background re- 
gion by ±5%, and derived two new foreground+CXB best-fit 
emission models. The best-fit parameters of these fits were 
then used to refit each ICM spectrum, while at the same 
time changing the normalization of the corresponding QPB 
spectra associated with the ICM bins by ±5%. MOS and 
pn QPB normalizations were raised and lowered simultane- 
ously. Figure [5] shows the changes on the best-fit parame- 
ters of each ICM bin, by comparing the new best-fit values 
with those summarized in Table IAT1 Again, the best-fit pa- 
rameters of bins 13 and 14 depend strongly on background 
uncertainties, and a 5% error in the normalization of the qui- 
escent particle background yields temperatures inconsistent 
with the best-fit results presented in Section 13.31 However, 
the best-fit parameters of the other bins are consistent with 
our previous results. 



4.3 Solar wind charge exchange 

Solar wind charge exchange (SWCX) photons, released in 
the aftermath of collisions between atoms and solar wind 
ions, sometimes affect X-ray observations below ~ 2 keV. 
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Figure 10. +7 /0 +6 ratio. Black points correspond to measure- 
ments taken during the time of the XMM-Newton observation of 
CIZA J2242. 8+5301, while the dashed and dotted lines show, re- 
spectively, the mean and maximum of these data points. 



The main lines appearing in the SWCX spectrum are Cvi 
(0.37, 0.46 keV), Nvi (0.42 keV), Ovn (0.56 keV), Ovm 
(0.65, 0.80 keV), Neix (0.92, 1.15 keV), and Mgxi (1.34 keV). 
Of these, the carbon and oxygen lines are generally the 
strongest (e.g. ISnowden et al.ll200j; ICarter fc Sembavll2008l ; 
ICarter et all l201lh . ISnowden et alj \2Q04 l has shown that 
periods of enhanced SWCX are associated with an increase 
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in the +7 /0 +6 ratio. In Figure 1101 we show the hourly 
+7 /0 +6 ratio over a period of 27 days (12 days prior and 
following the XMM-Newton observation). The solar wind 
data was obtained with the SWICS (Solar Wind Ion Com- 
position Spectrometer) instrument on board of the ACE 
(Advanced Composition Exporer) satellite. We only show 
FLAG=0 (i.e., good quality) data points. As Figure 1101 
shows, during the XMM-Newton observation the +7 /0 +6 
ratio was among the lowest recorded over almost a whole 
month. Therefore, SWCX emission is negligible. 



5 EVIDENCE FOR SHOCKS IN THE ICM 

In Figure [4l a surface brightness discontinuity can be seen 
to the south, at approximately 1.3 Mpc from the centre of 
the sectors in Figure [3] We modelled it assuming spherical 
symmetry, and an underlying density profile described by 
two power-laws with different normalizations and indices, 
one for each side of the discontinuity: 
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where C is the shock compression, r s h oc k is the cluster- 
centric distance at which the density discontinuity is ob- 
served, and n s hock is the electron number density of the 
post-shock gas. Starting from this density model, the surface 
brightness, Sx, was calculated by projecting the emissivity 
along the line of sight, 5*x = J^ L edZ, with the emissivity, e, 
defined as: 



e = n c mA(T, Z) 

where the ion density is rii 
function, 



(2) 

n e /1.2, and A is the cooling 



A(T, Z) = Ci(kT) a + C 2 {kTf + C 3 



(3) 



in units of 10 22 erg cm 3 s 
Ci = 8.6 x 10~ 3 ergcm 3 s~ 1 keV 17 
erg cm 3 s" 1 keV~° 5 , and C 3 



and with a = -1.7, /3 = 0.5, 
C 2 = 5.8 x 10" 2 
6.4 x 10~ 2 ergcm 3 s _1 for 
Z = 0.3Z Q (e.g., iRuszkowski et all 12004 and references 
therein). L is the maximum distance from the centre of the 
cluster, measured along the line of sight; since we assumed 
spherical symmetry, this distance is equal to the maximum 
projected distance to which the surface brightness profile 
is calculated. Because the number of points in the putative 
pre-shock region of the surface brightness profile is small, 
the power-law index of the density model on this side of 
the discontinuity was fixed to 2.466, which co rresponds to 
three times the average j3 index obtained by lEckert et al.l 
i|2012h when fitting a beta model to the profiles of 17 non- 
cool core (NCC) ROSAT clusters between 0.65 and 1.3r2oo- 
Furthermore, the distance of the jump from the centre was 
fixed to the average of the two radii on both sides of the 
observed discontinuity. In conclusion, the fits had three free 
parameters: C , n s h oc k, and 72. 

The southern surface brightness profile crosses the radio 
relic. However, assuming a shock is present at the relic, the 
expected pre-shock region (i.e., the region right in front of 
the relic) lies beyond 12 arcmin from the detectors' centre, 
and it is not included in our imaging analysis. Instead, the 
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Figure 11. EPIC spectra extracted from the (putative) post- 
shock region behind the northern relic, and the best-fit thermal 
model. MOS1 and MOS2 data are plotted in black and red, re- 
spectively, while pn is plotted in green. The bottom pannel shows 
the ratio of the data to the model. 



observed surface brightness discontinuity at ~ 1 Mpc is due 
to the fainter region SE of the relic, at angles roughly be- 
tween 240° and 270° . We assumed a temperature of 7 keV on 
both sides of the discontinuity, which is an acceptable esti- 
mate based on the temperature map in FigureJS] The best-fit 
to the surface brightness profile has \ 2 /A.o.i. = 13.9/19, and 
yields C = 1.35 ± 0.14, corresponding to a Mach number of 
1.24 ± 0.097. We tested the effect of varying 71 on the re- 
sultin g Mach number. The highest j3 derived bv lEckert et al.l 
|2012l ) is /3 = 0.963±0.054, for cool core (CC) clusters whose 
profiles were fitted with a beta model in the radial range 
0.65 - 1.3r 200 , while the lowest is P = 0.677 ± 0.002, for 
NCC cluster profiles fitted in the radial range — 0.7r2oo- 
Therefore, we varied 71 between 2.025 and 3.051, and found 
Mach numbers in the range of 1.2 — 1.3. If a Mach number 
of ~ 1.3 reflects the strength of the shock at the southern 
relic, it would be the lowest Mach number shock known at 
a relic. This would be particularly surprising so far from 
the merger centre, where the sound speed is expected to be 
low, especially with the southern relic having a counter-relic 
that is expected to trace one of the strongest merger shocks 
(M ~ 4 — 5 base d on the radio-der i ved M ach number for 
the northern relic, Ivan Weeren et all (|201Ch ). 

It is also interesting that the shock extends, to the east, 
beyond the souther n relic. A similar featu re has been pre- 
viously observed by iMazzotta et al.l (|201ll ) in RXJ1314.4- 
2515, in which a shock extends beyond the western relic to 
the north and south. It is possible that the shock strength 
varies across the shock boundary, and the Mach number of 
the putative shock across the southern relic is actually higher 
than 1.3. However, studying this would require further ob- 
servations, this time focused on the southern relic. 

In the north, the temperature peaks at w 9 keV near 
the western tip of the relic, as can be observed in Figure 
[U However, none of the bins provides an adequate repre- 
sentation of the (putative) post-shock gas. Therefore, we 
measured the temperature in a partial annulus extracted 
from the 53° — 123° sector shown in Figure [31 between radii 
of 0.75 and 1.35 Mpc; we note that this region overlaps 
most of the relic. The spectral model was the same as that 
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Table 4. Best-fit parameters to the "smudge" spectra. 





THERMAL MODEL 


NON-THERMAL MODEL 


r 

A/"x 


Q 1+7.7 
(1.4 ±0.12) x 10~ 5 


1.6 ±0.14 
3.6ir 6 8 x 10-" 


x 2 /d.o.f. 


86.88/96 


87.79/96 



described in Section 13.31 but the background parameters 
were kept fixed to the values listed in Table [3] The best- 
fit temperature was found to be 9.2i~ 2 'g keV, and the fit had 
X 2 /d.o.f. = 225.43/244. The fitted spectra are shown in Fig- 
ure llll Therefore, we see the temperature slightly increasing 
from 7 — 8 keV at the centre of the merger, to approxi- 
mately 9 keV behind the northern relic. The radio spectral 
index at the front of the relic predicts, in the linear regime, 
a Mach number of 4.6io'g. If this Mach number matched 
the density and temperature jumps across the relic, then we 
would expect a shock compression o f approximately 3.5, and 
a pre- shock temperature of ~ 1 keV. lAkamatsu fc Kawaharal 
|201ll ). using Suzaku observations of the cluster, measured 
a pre-shock temperature of approximately 1.6 keV. How- 
ever, their post-shock temperature is only 6.7 keV, much 
below the temperature measured behind the relic using our 
XMM-Newton observations. As also pointed out in Ogrean 
et al. (2012), most of this di fference is caused by an incorrec t 
background model used bv lAkamatsu fc Kawaharal (|201 ll ). 
who omitted the HF component. At this point, it is unclear 
how an additional foreground component would affect the 
Suzaku measurement of the pre-shock temperature. 

The low count statistics of our data do not allow us to 
constrain the density or temperature beyond the northern 
relic. However, the temperature map indicates an increase 
in temperature at the position of the relic. Such a tempera- 
ture profile towards merger shocks is not surprising. As the 
shock propagates from the centre towards the outskirts, the 
downstream gas furthest from the current shock position 
has had the longest time to cool. Assuming that the ICM 
can be approximated by an isothermal gas sphere, the ICM 
that was traversed m ost recently by the sh ock will have the 
highest temperature (|Roettiger et al.ll 19971 ). However, these 
approximations are simplistic, and substructure, projection 
effects, as well as the underlying cluster temperature profile 
can erase the temperature gradient. In fact, a positive tem- 
perature gradient in the direction of the relic has not been 
observed in the two relic-hosting mergers with confirmed 
shocks. 

In conclusion, we find evidence of a weak shock with 
M ~ 1.2 — 1.3 near the southern relic and, interestingly, the 
shock extends beyond the SE tip of the relic. In the north, 
the temperature is highest at the relic, peaking at roughly 
9 keV. Unfortunately, our XMM-Newton observations do not 
allow us to constrain the density or the temperature in front 
of the northern relic, and therefore we are not able to mea- 
sure the Mach number of the putative shock. 



Radio flux (mJy/beam) 




0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 



-i — | — , — , — , — | — , — , — , — | — , — , — , — | — , — , — , — r 

0.4 - 



0.2 - 



0.0 - - 
-0.2 - > 
-0.4 - 

Li i i i i i i i i i i i i i i i i i 

-0.4 -0.2 0.0 0.2 0.4 
Mpc 

Figure 12. GMRT 1280 MHz radio flux map, focused on 
the AGN at the eastern tip of the northern relic in CIZA 
J2242. 8+5301. While the AGN is positioned in continuation of 
the relic, it is detached from it, and the AGN jets are pointing 
towards NE, perpendicular on the relic. 

6 THE NORTHEASTERN AGN 

The XMM-Newton surface brightness map shows a faint, 
half-Mpc-long "smudge" located northeas t from the cluster. 
This c oincides with an AGN identified by Ivan Weeren et al.l 
l|201Cl ) in the radio. We extracted EPIC spectra in a par- 
tial annulus of area ~ 6 arcmin 2 overlapping the X-ray 
"smudge" , and fitted them with two different spectral mod- 
els: a non-thermal model described by a redshifted power- 
law, and a thermal model, described by an absorbed APEC 
model. The photon indices, temperatures, and normaliza- 
tions were coupled between detectors, while the background 
was fixed to the best-fit background model summarized in 
Table [3] The results of the fits are summarized in Table 3J 
and the fitted spectra are shown in Figure [13] The errors on 
the thermal model are quite large and tip the balance toward 
non-thermal emission from the smudge. However, we cannot 
clearly differentiate between thermal and non-thermal emis- 
sion. Unfortunately, the XMM-Newton data does not allow 
us to fit a combined thermal plus non-thermal model to the 
"smudge" spectra. 

It is not clear from the present observations whether 
the "smudge" is associated with the AGN or with the relic. 
In low frequency radio observations, the AGN appears de- 
tached from the relic, as can be seen in Figure[l2] Its jets are 
oriented perpendicular to the "smudge" , while the "smudge" 
is almost perfectly aligned with the southern boundary of 
the eastern tip of the northern relic. It is possible that a ra- 
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Figure 13. The "smudge" spectra, fitted with a thermal, and non-thermal model. MOS1 and MOS2 data are plotted in black and red, 
respectively, while pn is plotted in green. The bottom pannel shows the ratio of the data to the model. 



dio galaxy is ram-pressure stripped as it is moving through 
the ICM, possibly driven by the motions triggered by the 
merger, and this might be what we are observing in the X- 
ray. However, if this is the case, it is difficult to explain why 
at least one jet is so clearly oriented perpendicular to the 
apparent direction of motion. 

If the "smudge" was caused by a merger shock, we 
would be more likely to see it at the front of the relic, and 
it is unclear why the emission would be so enhanced only 
on one side of the shock. One possibility is the existence of 
a cosmic filament oriented roughly NE-SW. The density in- 
side this filament would be higher than the density of its sur- 
roundings, which would increase bremsstrahlung emission in 
the part of the shock that crosses the filament. The shock 
speed along the filament would also be lower, which would 
curve the eastern tip of the relic inwards, and could diffuse 
the synchrotron-emitting particles; indeed, the fact that the 
geometrical centre of the relic system, assuming sphericity, 
is shifted slightly to the west, and that the radio structure 
S-SE of the relic is more filamentary could be an indication 
of this geometry. On the other hand, the magnetic field in 
the ICM region crossed by the filament would also be higher 
than in its surroundings. The synchrotron power scales as 
P S yn oc B 2 , where B is the magnetic field. The radio surface 
brightness of the northern relic in CIZA J2242. 8+5301 is 
surprisingly uniform, therefore, assuming a filament crosses 
the eastern part of the relic, the lower shock speed inside the 
filament should compensate for the higher magnetic field. 



7 THE MERGER SCENARIO 



Numerical simulations by Ivan Weeren et al.l ()201lh predict 
that CIZA J2242. 8+5301 is a merger occurring in the plane 
of the sky, between two almost equal-mass clusters (mass 
ratio 2 : 1). This simple geometry minimizes complexities 
introduced by projection effects, and would indicate that the 
cluster is a textbook example for studying cluster mergers 
and large-scale shocks. However, also from numerical sim- 
ulations, we know that such a merger scenario would re- 
sult in two outward-moving shocks (advancing through the 
ICM perpendicular to the merger axis), two radio relics trac- 



ing these shocks, and a relatively symmetric temperature 
and pressure distribution. While the pressure map of CIZA 
J2242. 8+5301, shown in Figure [6] is roughly the same on 
both sides of the merger axis, the temperature map is far 
from symmetric. Two hot bins enveloped in cooler plasma 
are seen north and south of the centre. The existence of 
these hot spots is indeed predicted by numerical simula- 
tions of cluster mergers, which show that the cores of the 
merging clusters can be heated to temperatures above the 
post-shock temperature, during the late phases of the merger 
event (e.g.. Ivan Weeren et al.|[201ll ). However, to the east, 
the temperature increases from the centre into a hot "wall" 
with temperatures of 9 — 12 keV. The existence of this hot 
region is puzzling, as the gas is not associated with a radio 
relic, and we do not identify a surface brightness discontinu- 
ity at its position. Furthermore, the presence of a hot eastern 
"wall" is in apparent contradiction with the simple merger 
geometry indicated by numerical simulations, suggesting a 
more complicated merging scenario, or a lack of understand- 
ing on our part of the complex structures formed during real 
cluster mergers. 

One possibility is that the merger event had a small 
impact parameter. Th is scenario resemb l es tha t shown in 
fig. 2 (t » 2.4 Gyr) of Ivan Weeren et alj |201ll ). if a small 
counterclockwise rotation is applied to the image, such that 
the relics are aligned along the N-S direction. In this merger 
geometry, however, a "wall" of cold gas should be present 
east of the merger a xis, roughly 0.5 Mpc from the merger 
centre (fig. 3, right, of lvan Weeren et al.l (|201ll 'l). This is the 
opposite to what is observed in the temperature map, as 
the gas 0.5 Mpc east of the cluster centre is in fact very hot. 
If the "wall" of cold gas predicted by numerical simulations 
actually forms in real mergers, then in seems more likely that 
it is present west of the centre of CIZA J2242. 8+5301, while 
to the east, the ICM is heated by a process or event that is 
unaccounted for in the numerical simulations (e.g., a triple 
cluster merger, with two clusters colliding in the plane of the 
sky, and a third one moving roughly along the line of sight). 
If this is the case, then the main merger event involved a 
collision between two almost equal mass clusters (mass ratio 
between 1:1 and 2:1), with the less massive cluster entering 
the gravitational well of the more massive cluster from the 
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SE. There are indications that this scenario is also supported 
by the surface brightness distribution, which appears slighly 
broadened in the SW-NE direction, compared to the SE- 
NW. The constraints set by the X-ray morphology on the 
merger scenario will be discussed in more detail in Ogrean et 
al. (in prep.), using high resolution maps from a deep (200 
ks) Chandra observation. 

To summarize, although we are not able to measure 
the temperature beyond the spectacular northern relic, the 
temperature distribution to the north and south of the clus- 
ter centre is strikingly similar to that predicted by numerical 
simulations of merging clusters triggering large-scale shocks. 
However, while for CIZA J2242. 8+5301 numerical simula- 
tions suggest a "simple" merger between two nearly-equal 
mass clusters colliding in the plane of the sky, this scenario 
is not supported by our X-ray results. We observe a very 
asymmetric temperature distribution, with a "wall" of hot 
plasma separating the merger centre from the eastern-most 
outskirt regions. The most likely merger scenario is that in 
which a less massive cluster entered the gravitational well of 
the more massive cluster from the SW, with a small impact 
parameter. 



8 SUMMARY 

Based on radio data and numerical simulations, CIZA 
J2242. 8+5301 appears to be a textbook example of a merg- 
ing galaxy cluster hosting a double-relic system. The pre- 
dicted merger geometry is simple: two clusters of almost 
equal masses merging along the north-south direction and 
essentially in the plane of the sky, which collid ed with a 
small impact parameter (|van Weeren et al1l201lh ■ However, 
this simple picture is only partially supported by the XMM- 
Newton results presented above. Below we summarize our 
main findings: 

• In the X-ray, the cluster does indeed show an extremely 
disturbed, N-S-elongated morphology, indicative of a merger 
along the N-S direction. The southern part of the main X-ray 
structure is linear and brightest toward its southern tip, rem- 
iniscent of a bullet-like morphology. From the N-NE tip of 
this linear structure, the emission curves in a west-pointing 
arc of enhanced surface brightness. This X-ray morphology 
is more complex than predicted by numerical simulations. 

• Because of the relatively high instrumental background 
level of XMM-Newton, we were not able to characterize the 
surface brightness profile beyond the spectacular northern 
relic, and therefore find no evidence of shock compression 
near the relic. To the south, however, slightly east of the 
southern relic, we find evidence of a shock of Mach number 
~ 1.2-1.3. Due to the size of the XMM-Newton FOV, we do 
not have data immediately beyond the relic. Nevertheless, 
we speculate that a shock is present at the southern relic, 
and that what we observe is its extension beyond the de- 
tected boundaries of the relic. Possible reasons for different 
relic and shock lengths include Mach number fluctuations 
across the shock, coupled with weak cosmic ray injection 
efficiency at low Mach number shocks, or simply detector 
limitations. 

• While numerical simulations of mergers with low mass 
ratios and low impact parameters predict relatively symmet- 
ric temperature distributions, this is not what is observed 



in CIZA J2242. 8+5301. The temperature increases at the 
northern relic, as expected for shocks at which the increase 
in gas temperature dominates over the underlying temper- 
ature profile of the cluster. A smaller increase in temper- 
ature is also visible near the southern relic. What is more 
surprising is a "wall" of hot gas located east of the clus- 
ter centre. This "wall" is not associated with a radio relic. 
However, it appears to extend into the hot region behind the 
northern relic, forming approximately a semi-ellipse of hot 
plasma around the cluster. This temperature distribution is 
not what numerical simulations predict for a simple merger 
geometry. 

• A "smudge" of enhanced X-ray emission can be seen 
aligned with the easter tip of the northern radio relic. Part 
of the smudge also overlaps the position of an AGN visible in 
the radio, whose jets point away from the cluster centre, per- 
pendicular to the relic. It is not clear whether the "smudge" 
is associated with density compression at the relic, or with 
the AGN. Its spectrum appears to be better described by 
a non-thermal model. However, given the orientation of the 
AGN jets and the fact that it is clearly detached from the 
relic, we consider that the "smudge" is more likely a result 
of shock compression at the relic. The fact that it is visible 
only along ~ 1/5 of the length of the relic could be the result 
of cosmic structure (e.g., a cosmic filament), or a sign that 
the properties of the (putative) shock at the northern relic 
vary strongly across the relic's length. 
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Figure Al. Bin map used for creating the temperature and pres- 
sure maps in Figures [6] and IB 1 1 Best-fit spectral parameters for 
each bin are listed in Table |A"T1 and the fitted spectra arc shown 
in Figures lAll and lB2l 



Table Bl. Fitted background and foreground parameters, ob- 
tained using a fixed X-ray column density of Nn = 3.22 X 10 21 
cm . 
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APPENDIX A: SPECTRAL FITS 
APPENDIX B: RESULTS WITH FIXED iV H 
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Table Al. Best-fit parameters to the spectra extracted from the bins in Figure [All Temperatures are in units of keV, normalizations 
are in default Xspeo units per square arcmin, and X-ray column densities are in units of 10 21 cm -2 . The bins were fitted in parallel. 
The fit with a fixed JVh of 3.22 X 10 21 cm -2 had a reduced x 2 of 1.14. The colums show: (1) bin number, (2) best-fit temperature for 
iV H = 3.22 X 10 21 cm -2 , (3) best-fit normalization for JV H = 3.22 X 10 21 cm -2 , (4) best-fit temperature obtained with a free JV H , (5) 
best-fit normalization obtained with a free Nn, (6) best-fit Nn. 
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Figure Bl. Temperature and pseudo-pressure maps of the cluster, obtained using a fixed X-ray column density of TVjj = 0.322 X 10 21 
cm -2 . Spectra were extracted from bins that follow the surface brightness contours, and have a total of ~ 3600 counts after the 
subtraction of the instrumental background. The pressure was calculated as the product between the temperature and the square root 
of the normalization of the spectral component describing 1CM emission. Overlaid are the radio contours also shown in Figure [T] 
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